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Abstract. We describe the horizon of a quantum black hole 
in terms of a dynamical surface which defines the boundary of 
space-time as seen by external static observers, and we define a 
Ph! path integral in the presence of this dynamical boundary. Using 

J|h I renormalization group arguments, we find that the dynamics of 

the horizon is governed by the action of the relativistic bosonic 



membrane. From the thermodynamical properties of this bosonic 
membrane we derive the entropy and the temperature of black 
holes, and we find agreement with the standard results. With 
this formalism we can also discuss the corrections to the Hawking 
temperature when the mass M of the black hole approaches the 
Planck mass Mpi. When M becomes as low as (10 — 100)Mpj 
a phase transition takes place and the specific heat of the black 
hole becomes positive. 



In a recent paper [ED] we have put forward a proposal for a quantum 
description of black holes. Following the work of 't Hooft and of Susskind 
and coworkers ||, we have considered the black hole horizon, as seen by 
an external static observer, as a dynamical surface endowed with physical 
degrees of freedom, and we have suggested that the dynamics of this surface is 
governed by an action principle. The simplest action describing a relativistic 
bosonic membrane in 3+1 dimensions is [|J 

-rf^V=h, (i) 



-'memb 

where T is the membrane tension, and = (r, er^o^) parametrizes the 
membrane world-volume; h is the determinant of the induced metric hij = 
g fJlU diC' J 'dj( l/ and x^ = C M (£) gives the embedding of the membrane in 3+1 
dimensional space-time. The target space metric g^ u is taken to be the 
Schwarzschild, Reissner-Nordstrom, or Rindler metric. We limit our con- 
siderations to non-rotating black holes. 

Classically, the action (1) describes a membrane which moves in the black 
hole background, approaching the horizon asymptotically. For instance, in 
the case of the Rindler metric ds 2 = —g 2 z 2 dt 2 + dx 2 + dy 2 + dz 2 we have 
found a classical solution of the form 

*(r) = -£-, (2) 
cosh gr 

which approaches asymptotically the classical horizon, z = 0. We have dis- 
cussed the quantization of the action (1) in a minisuperspace approximation, 
and we have found that, for a given value of mass and charge, the quantum 
state of a black hole is not uniquely determined, but rather there exists a 
quasi- continuum of levels corresponding to excitations of the membrane. In 
the classical description the state of a (non-rotating) black hole is charac- 
terized only by the "macroscopic" parameters M, Q, and therefore a coarse 
graining over the membrane levels is implicit. The membrane approach there- 
fore provides a microscopic explanation of the black hole entropy. This de- 
scription also suggests a quantization of the area of the horizon, in agreement 
with various heuristic arguments existing in the literature ||. 

In this Letter we examine the membrane description from the path inte- 
gral point of view. We discuss how the action (1) emerges from a renormal- 
ization group analysis, and we show how the thermodynamical properties of 
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black holes can be computed in our approach. We will find that the black 
hole radiation can be understood in terms of fluctuations of the horizon - a 
point of view stressed in particular by York 0. 

Let us consider the definition of the path integral in the presence of black 
holes. From the point of view of a fiducial observer^ only the degrees of 
freedom outside the horizon can influence the dynamics. So, in a path inte- 
gral approach, we would like to perform the integration only over the field 
variables g^x) with x outside the (apparent) horizon. However, a com- 
plication arises immediately, since the position of the horizon is determined 
by the metric itself: when the metric fluctuates, the horizon fluctuates (see 
e.g. H D, U), and as a consequence the number of variables g^(x) which we 
would like to include as integration variables in the path integral changes, 
and we are faced with the problem of defining a path integral in which the 
number of integration variables fluctuates. A possibility would be to proceed 
to a gran canonical ensemble.0 However, we rather proceed as follows. First, 
we introduce a fixed spherical surface located at r = r + + e, where r + is 
the radius of the classical horizon and e is larger than the typical quantum 
fluctuations of the horizon, which are on the order of a few Planck lengths. If 
r > r + +e the corresponding variables g^ v are inserted as integration variables 
in the path integral. The question is what to do in the shell r + < r < r + + e. 

In general, when evaluating a path integral in field theory, a crucial point 
is the identification of the physically relevant variables, as opposed to fast 
varying variables which can be integrated out, in the spirit of Wilson's renor- 



malization group [10]. In a non-perturbative regime, the physically relevant 
variables are very different from the fields which appear in the Lagrangian. 
An obvious example is provided by QCD at length scales on the order of 
one fm, where the physical degrees of freedom are hadrons rather than free 
quarks and gluons, and the gluon field between a quark-antiquark pair is 
squeezed in a flux tube whose dynamics is governed, at the effective level, by 
the action of a non-critical relativistic string. In our case, a similar situation 
arises because, from the point of view of a fiducial observer, the region within 
a few Planck lengths from the classical horizon is a region of super-Planckian 



1 For a non-rotating black hole a fiducial observer is defined as a static observer outside 
the horizon. The precise definition of fiducial observers for rotating black holes, together 
with a detailed description of the membrane paradigm for classical black holes, can be 
found in ref. 

2 I thank Enore Guadagnini for this remark. 
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temperatures and very large fields. 

It is natural to identify the "collective variables" which play the main 
role with the coordinates (^(6, <p, t) which define the position of the horizon, 
x n _ £m Qf course; the variables C 1 are uniquely fixed by the metric g^ v . 
So, the integration variables g^ u (x), restricted by the condition that x lies 
between the horizon, as determined by g^ u itself, and the surface r = r + + e, 
are decomposed into the physically relevant variables £ M plus "fast variables" . 
The effective action for the £ M is obtained integrating over the fast variables. 
While such integration would be very hard to perform explicitly, the general 
form of the effective action is fixed by invariance principles. The simplest 
term which can arise is just the membrane action, eq. (1), with a membrane 
tension T which is in principle derivable from the underlying theory (even 
if its actual calculation is a difficult non-perturbative problem, analogous to 
the computation of the string tension in QCD). 

Thus, we can write the partition function as 

V I Tin J 1 



V(» Vg„„e*\ (3) 

where Q is the region of space-time outside the sphere r = r + + e. The action 
I is given by 

I ^memb -^grav i (4) 

where / mem b is the membrane action, eq. (1) (plus higher-order terms which 
will be in general generated by the renormalization group procedure), and 
igrav is the gravitational action of the variables g^x) with x G Q. A crucial 
point is that, since the region Q has a boundary dQ, we must consider the 



gravitational action supplemented with the boundary term [11, 12 



I = / d^x^gR-^-l dPxy/^K, (5) 
16n Jn 8n Jan 

where K is the trace of the extrinsic curvature of the boundary dQ, i.e. of 
the surface r = r + + e, and 7 is the determinant of the induced metric on 
this surface (not to be confused with h which denotes the determinant of the 
induced metric on the dynamical surface x^ = to avoid confusions let us 
also stress again that dQ is a fixed mathematical surface while x^ = C M (£) 
gives the position of a physical, fluctuating, surface). 

As it stands, eq. (^|) is only formal, since we have not specified what we 
mean by integration over g^. The problem of the definition of the path 
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integral over g^ has been first discussed in the well-known papers by Gib- 



bons and Hawking |12| and Hawking [13]. However, the important quantum 



fluctuations are that in the region r + < r < r + + e, which have already been 
taken into account by the renormalization group procedure. Thus, the dif- 
ficult problems connected with a proper definition of the path integral over 
the metrics can be avoided: to lowest order in the semiclassical expansion, 
the variables g^ u at r > r + + e can be replaced by the classical metric g^ u . 

The definition of the path integral is completed performing a Wick rota- 
tion on the membrane world-sheet, r — > it. (Note that r appears both in d 4 x 
and in d 3 x in eq. (H) so that both terms pick a factor i, and 7 becomes the 
Euclidean induced metric.) Performing the Wick rotation only at this stage, 
we avoid the problems connected with the definition of Euclidean quantum 
gravity. Thus, at this order, we write 

e- J — /pC M e- J — h )\ g ^= g « ■ (6) 

The path integral in the above equation is the partition function of an Eu- 
clidean membrane theory. Note that in the classical background metric g c ^ v 
the curvature outside the horizon is zero and only the extrinsic curvature 
term contributes to / gra v 

In order to compute the thermodynamic properties of black holes in the 
membrane formalism we proceed as follows. First, it is convenient to choose 
the gauge (° = r in the action (1); thus, for a black hole metric r coincides 
with the (Euclidean) time of asymptotic observers. The solution given in 
eq. (0) refers to this choice of gauge. Then we consider the partition function 
over fields periodic in r with period /3oo and we evaluate the partition function 
in saddle point. In the case of the Rindler metric a solution of the classical 
equations of motion is given by eq. (§), and in Euclidean space it becomes 
an instanton of the form 

z(t) = -*»-. (7) 
cos gr 

We see that this solution is periodic in r with period (3^ = 2n/g. So it 
contributes to the path integral for the partition function of a canonical 
ensemble with temperature at infinity 

r-£, (8) 
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which is the well known result for the temperature measured in Rindler 
space 
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Let us now discuss the Schwarzschild black hole. The equation of motion 
of the membrane before the rotation r — > ir is Iffll 



ra' 



rr + 2(a 2 - r 2 ) + — (a 2 - 3r 2 ) = 0, (9) 
2a 

where a = 1 — 2M/r and a' = da /dr. In more geometrical terms the above 
equation can be written K [Q = 0, where K [Q is the trace of the extrinsic 
curvature of the surface x^ = C M (0- We define 

z = 4M (1 - — f/ 2 . (10) 
r 

This definition is motivated by the fact that, with this change of variables, 
together with x = 2M(6 — (7r/2)) and y = 2M<p, the Schwarzschild metric for 
a black hole with large mass goes into the Rindler metric, plus corrections 
of order 0(z 2 /M 2 ),0(x 2 /M 2 ), with the constant g which appears in the 
Rindler metric identified with the surface gravity k = 1/(4M). For the sake 
of comparison with the Rindler limit, in the following we will denote the 
surface gravity by g. Note that < z < 1/g as 2M < r < oo. In terms of z 
the equation of motion (|Sp reads 

zi - 2z 2 + g 2 z 2 {\ + 3#V)(1 - g 2 z 2 f = 0. (11) 

li g ^ this equation reduces to the equation of motion of the membrane in 
Rindler space, z'z — 2z 2 + g 2 z 2 = 0, which has the solution given in eq. (0). 
This means that, if the mass of the black hole M goes to infinity, then 
g = 1/ (AM) — > and equation (|TTJ) (after the Wick rotation r — > ir) admits 
a periodic solution with period 2n/g. Thus, for black holes in the limit 
M — > oo, we recover the Hawking temperature ||15|], 



T=-^— . (12) 

Beside reproducing the Hawking temperature in the large mass limit, eq. (|i"T|) 
allows us to study how the result is modified at finite mass. The semiclassical 
approach used in the standard computation of black hole radiation is valid 
only if M 3> Mp\, where Mpi is the Planck mass, set equal to one in our units. 
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Thus, in general, finite mass corrections should be expected. The situation, 



however, is more complicated. Eq. (11) cannot be studied expanding pertur- 



batively z(r) = z (r) + g 2 z 2 (r) + g A z i {r) + . . . because z (t) = a/ cos(gr) 
diverges at gr = tt/2 and therefore the non-linear terms in eq. (|ll|) are not 
small, even if g <C 1 (unless g is strictly zero). 

To study eq. flllD we introduce the function w(r) = gz{r) and the variable 
t = gr. In Euclidean space, the equation for w(t) is 

W 1F ~ 2( ^ )2 ~ ^ 2(1 + ~ ^ f = ° • (13) 

By definition, < w < 1 since < z < 1/g. Apparently, all dependence on g 
has disappeared. However, when solving eq. flT3| ) with the initial conditions 
u>(0) = a, w'(0) = 0, where a is a constant, < a < 1, we will find peri- 
odic solutions, see below. The period of these solution in general depends 
on the amplitude a, because of the non-linearity of the equation. Thus, one 
should ask what fixes the typical scale of a. In the path integral, a appears 
as a collective coordinate of the instanton solution. In principle, the distri- 
bution of values of a is fixed by the full theory, similarly to what happens 
for the instanton radius in gauge theories. Qualitatively, it is clear that the 
typical fluctuations of the horizon are of order one in Planck units, so the 
main contribution comes from solutions with z(0) ~ 1 and therefore a ~ g. 
Fluctuations of the horizon with larger values of a should be suppressed, and 
anyhow cannot be computed in the membrane formalism because, as we have 
discussed in ref. fluctuations with amplitude a ^> 1 correspond to highly 
excited state of the membrane, for which our description is not self-consistent 
since we are neglecting the backreaction of the membrane on the metric or, 
equivalently, we are neglecting higher order terms in the membrane effective 
action. 

Thus, the dependence on a that we will find below can be translated, apart 
from numerical constants, into a dependence on g just identifying a ~ g. In 
this way it is possible to discuss, at least qualitatively, finite mass corrections 
to the Hawking temperature. 

We have integrated numerically eq. (|1^) for different values of a. Let 
us introduce the function u(t) = l/w{t). The results are as follows. For 
very small a (a < 0.004) we find a periodic solution u{t) which is essentially 
undistinguishable, within our numerical accuracy, from the Rindler limit with 
the same initial conditions, i.e. from the function (1/a) cost. At a = 0.004 
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we can put an upper bound on the difference between the period P of this 
solution and the period Pq = 2n of the unperturbed solution, AP/P < 1CT 3 . 

Identifying the period of the solution with the inverse temperature of a 
black hole with mass M such that g = 1/(4M) ~ a, we see that, even for 
black holes with mass as low as I /(AM) ~ 0.004, or M ~ (10 - 100)M Pb the 
corrections to the Hawking temperature are numerically neglegible - if they 
exist at all. 

For larger values of a a new phenomenon takes place. An instability 
occurs, and the solution, instead of a form qualitatively similar to (1/a) cost, 
takes a form qualitatively similar to (1/a) | cost|, i.e., the negative part of the 
solution is flipped to positive values. For a close to this critical value a c ~ 
0.004, the numerical integration becomes unreliable, with some half periods 
flipped and others not flipped, randomly. The reality of the phenomenon is 
however clearly seen increasing a further. Then, in correspondence with the 
point t ~ 7r/2 at which the derivative changes sign, the value of uit) moves 
away from zero and instead of a cusp the function u(t) has zero derivative in 
t = i, and the numerical integration becomes stable again. The solution has 
now the form shown in fig. 1. The important new aspect is that this solution 
has a period which is one half as before, corresponding to a temperature 
which is twice as large as the Hawking temperature. 

This result shows the existence of a non-analyticity in g. A conservative 
interpretation of this non-analyticity is that it fixes the limits of validity of 
the membrane description, or more in general of the semiclassical approach. 
Certainly this is not surprising since we know that the semiclassical expansion 
must break down when M is not large compared to one. We find that this 
happens when 1/(4M) ~ 0.004, or M ~ 10 — 100, a reasonable result. 

Of course, it is very tempting to see what our equation predicts for even 
larger values of a, corresponding to an even smaller value of the mass - with 
the remark that it is by no means obvious that our results can be estrapo- 
lated beyond this point of non-analyticity. We see that the equation suggests 
that, at the critical point, the black hole still behaves as an approximately 
thermal object, with a temperature twice as large as the Hawking temper- 
ature. Increasing a further, the period increases smoothly (see fig. 2) and 
goes to infinity as a — > 1, corresponding to a black hole temperature which, 
after jumping from l/(8irM) to l/(47rM), starts to decrease smoothly, until 
it becomes zero when the mass is on the order of the Planck mass. A numer- 
ical fit for the period P as a function of a shows that, when a — > 1, the data 
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are very well reproduced by P = const x (1 — a) -1 ; identifying again a ~ eg, 
with c ~ 1, we get an expression for the temperature when M becomes on 
the order of M P j, 

T = Ag(l - eg) , (14) 
where numerically we find A ~ 1.5. The temperature formally vanishes for 

g = 1/c ~ 1. 

Independently of the numerical details, the important new qualitative 
feature is that, for masses below a critical value, the specific heat becomes 
positive. At the critical point, a second order phase transition has taken 
place. 

Although quite intriguing, the scenario for M below the critical value at 
present is only speculative and we do not discuss it further here. Let us come 
back to the safer domain of large black hole masses and discuss the black hole 
entropy. 

Having fixed (3^ at the standard value, we can now compute the entropy 
of black holes, evaluating the partition function Z, and therefore the free 
energy. The leading term is given by the extrinsic curvature term evaluated 
at the boundary dQ. At this point the computation is formally identical to a 



computation already performed by York ||16|| : the trace of the extrinsic curva- 
ture of a static spherical surface with a generic radius r in the Schwarzschild 
background is 

i^-V^-M*- 1 / 2 , (15) 



The square root of the determinant of the induced metric is = r 2 a 1 ^ 2 sin 6 
so that 

J grav = 12nM 2 - 8nMr . (16) 
The local inverse temperature (3 is defined as 

f3 = f3 0O a l l 2 = MM(l- 1 -^ . (17) 
Considering J grav as a function of r and (3, with M = M(r, (3) defined by 



eq. (|T7|), one finds |TJ| the entropy 

~d[3 



S = P(^]-I^ = ^M 2 , (18) 
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which agrees with the Bekenstein-Hawking value. It should be emphasized 
that this result is independent of r, the radius of the surface on which the 
extrinsic curvature is computed. This explains the agreement between the 
value of the entropy in the membrane approach and the result of the com- 
putation in Euclidean quantum gravity (12], [13], [16). In the latter case the 
boundary term is computed on a surface at infinity, while in our approach it 
is computed at r = r + + e. However, since the result is independent of r, the 
two methods give the same answer. The fact that the result is independent 
of r has also the important consequence that the dependence on the exact 
position of the surface dfl disappears. 

It is interesting to compare this derivation of the black hole entropy with 
a recent result of Frolov and Novikov [IT]]. There is an important conceptual 
difference, since in ref. ]17) the entropy is attributed to the modes of physical 
fields which are located inside the horizon and therefore invisible for a distant 
observer. Instead in our derivation, in the spirit of the principle of black 
hole complementarity ||, we have avoided any reference to the "inside" of 
the black hole. However, we think that some similarities exist between the 



approach presented in this paper and the approach of ref. [[L7j i i n both cases 
the entropy is obtained tracing over some degrees of freedom. In our case, 
the extrinsic curvature term which is responsible for the entropy is generated 
by the renormalization group procedure, which is a way of tracing over short 



scale fluctuations in the vicinity of the horizon, while in JTTj a distinction 
between "visible" and invisible" modes is made, and the entropy is obtained 
performing the appropriate trace over invisible modes. It turns out that the 
main contribution to the entropy is given by the "invisible" modes which 
propagate in the narrow shell r + < r < r + + e, with e on the order of a few 
Planck length. So, if we rather define the "invisible" modes as the modes 
which propagate in this shell, without making any reference to the inside of 
the black hole, the result S ~ A found in JljJ still goes through, but the 
derivation is now conceptually very similar to what we have presented. The 
proportionality constant between the entropy S and the area A can only be 
estimated, but not exactly computed, in the approach of ref. []T7J. In the 



membrane approach, instead, it can be computed exactly, as we have seen. 
The difference is due to the fact that we perform the "tracing" at the level 
of the action, via the renormalization group. This produces the extrinsic 
curvature term on the surface r = r + + e, with a coefficient which is uniquely 
fixed by the fact that it must cancel the term with second derivatives in the 



9 



Si 



inst 



Einstein-Hilbert action. 

We can also try to go further an compute the corrections to the lead- 
ing term in the entropy, inserting the instanton solution into the membrane 
action (and integrating over the collective coordinates of the instanton). Fur- 
ther corrections are obtained computing the fluctuations around the instan- 
ton. However, in general this produces divergent terms. For instance in 
Rindler space, where the instanton has the form ((?]), the action of the solu- 
tion is divergent, 

r 2ir I 

T'z dt—- (19) 
JO cos 2 1 

where T = T J dxdy is proportional to the area of the membrane and the 
proportionality constant, given by the integral in dt, diverges. Thus, in 
general we have (writing explicitly also the Newton constant G) 

S=±(1 + C)A (20) 

where A is the area and C a formally divergent constant. This divergence 
does not come out as a surprise. A similar divergence has been found by 
't Hooft p[ computing the number of states of a scalar field in the vicinity 
of the horizon. Its relation with the ultraviolet problem in quantum gravity 
has been discussed in where it has been suggested that this divergence is 



reabsorbed by the renormalization of Newton's constant. In the membrane 
approach, this divergence can be traced to the fact that, in the semiclassical 
approximation, the quantum membrane dynamics is governed by a potential 
which admits a continuum spectrum [I]]. However, this only happens if we 
extrapolate the membrane effective potential to the region very close to the 
horizon, where from the point of view of the fiducial observer the tempera- 
ture becomes super-Planckian and quantum gravity effects set in. Thus, the 
divergence in the entropy is actually a short distance problem, related to the 
fluctuations of the metric at the Planck scale Lp\. Instead, the membrane ap- 
proach that we have presented, being derived from a renormalization group 
procedure, can only give an effective description valid at distances larger than 
Lpi. 
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